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Abstract 

X-ray free-electron lasers (XFELs) may allow to employ the single particle imaging (SPI) method 
to determine the strnctnre of macromolecnles that do not form stable crystals. Ultrashort pnlses of 
10 fs and less allow to outrun complete disintegration by Coulomb explosion and minimize radiation 
damage due to nuclear motion, but electronic damage is still present. The major contribution to 
the electronic damage comes from the plasma generated in the sample that is strongly dependent 
on the amount of Auger ionization. Since the Auger process has a characteristic time scale on 
the order of femtoseconds, one may expect that its contribution will be significantly reduced for 
attosecond pulses. Here, we study the effect of electronic damage on the SPI at pulse durations from 
0.1 fs to 10 fs and in a large range of XFEL fluences to determine optimal conditions for imaging 
of biological samples. We analyzed the contribution of different electronic excitation processes and 
found that at fluences higher than lO^^-lO^® photons///m^ (depending on the photon energy and 
pulse duration) the diffracted signal saturates and does not increase further. A significant gain in 
the signal is obtained by reducing the pulse duration from 10 fs to 1 fs. Pulses below 1 fs duration 
do not give a significant gain in the scattering signal in comparison with 1 fs pulses. We also study 
the limits imposed on SPI by Compton scattering. 

PACS numbers: 61.80.Cb,61.80.Az,87.53.Ay,87.64.Bx 
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I. INTRODUCTION 


Modern x-ray crystallogr^ 
crystals at atomic resolution [l, 


w methods make it possible to determine the structure of 
2| . In the age of molecular biology, one of the most important 
questions in life science is the determination of protein structures. Unfortunately, many 


protein macromolecules, e^ecially membrane proteins, either do not crystallize or form only 
extremely small crystals |3|. Imaging such nanocrystals, or even individual macromolecules, 
is impossible with conventional x-ray sources, since the sample is destroyed before a high 
resolution diffraction pattern can be obtained jd]. 

These difficulties can be circumvented with the use of x-ray free-electron lasers (XFELs) 
7|. The ultrashort, intense, and coherent XFEL pulses can be used to obtain the diffrac¬ 
tion pattern of a small, sub-micron sample before it is destroyed by a Coulomb explosion j^. 
Conventional radiation dose limits jd] are surpassed by orders of magnitude, since the dom¬ 
inant damage mechanism, the breaking of chemical bonds, can be outrun by short XFEL 
pulses. In addition, for reproducible samples, many single shot patterns with random sample 
orientations can be taken and oriented to produce a full three-dimensional (3D) diffraction 


pattern 


12| . Using phase retrieval algorithms 


0 , 0 , 


a 3D image of the samples elec¬ 


tron density can be reconstructed. This approach is conventionally called nowadays single 
particle imaging (SPI) jl5| . 

In spite of signihcant progress in imaging biological samples at XFELs 16l-ll9|. experi¬ 
ments have highlighted severe challenges for such single-particle imaging experiments. Par¬ 


ticles tend to show high conformational heterogeneity 20, l2l|], and the coherent scattering 


signal to background ratio is 
even more powerful pulses 


m- 


ow. To improve the signal level it has been suggested to use 


However, the intensity of the x-ray pulse cannot be increased without limit. At a cer¬ 
tain threshold, extremely intense x-ray pulses strip all electrons from the atoms, leaving a 
sample that does not scatter any more. In addition, statistical fluctuations in the random 
ionization of the atoms due to the quantum mechanical nature of these processes produce 
a background signal that dominates the diffraction pattern for strong ionization. In our 


previous publication 


22| . we have shown that x-ray induced electronic damage limits the 


pulse fluence that can be reasonably employed in an experiment. 
It was also realized 


22 


23| that a major ionization mechanism in SPI experiments is 
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a secondary ionization by trapped Auger electrons. At the same time, Auger decay has 
a lifetime of several femtoseconds, which points to the possibility to outrun the impact 


ionization by trapped electrons with extremely short, possibly attoseconc 


was recently suggested that such pulses can be produced in principle 


m 


XFEL pulses. It 


25|. 


The aim of this paper is twofold. On one hand, we study to what extend the use of 
ultrashort XFEL pulses reduces the electronic damage of a typical biological sample. In 
particular we are interested in pulse durations from 0.1 fs to 10 fs, since this is the range 
where the suppression of the Auger decay is expected. On the other hand, we analyse 


an extension of the model used in 22|] by including additional ionization processes such as 
shake-off and Compton scattering. We also calculate the contribution of Compton scattering 
to elastically scattered radiation measured on the detector (see also recent work 2g where 
Compton scattering from a carbon cluster was analyzed). 


II. THEORY 


A. Elastic scattering 


In the following we recall the description of a single particle coherent diffraction exper¬ 
iment in the frame of kinematical approximation (see for details Ref. 22| and Appendix 

In a typical experiment a large number of single-shot diffraction patterns at different 
orientation of the particles will be measured. In the following we will assume that all 
these diffraction patterns can be perfectly aligned and averaged. If the incident beam is 
fully coherent over the sample area and has a uniform intensity distribution, we obtain the 
following expression for the averaged scattered intensity 221 

^ r 

< ^(q) >= / dtJ{t) < /*(q,t)/fc(q,t) > . (1) 


i,k=l 


Here, q is the scattering vector, J{t) is the intensity of the incoming pulse, and /i(q, t), R* 
are the time-dependent form factor and position vector of the Ath atom, respectively. The 
brackets (...) denote averaging over many pulses, and the summation is performed over all 
atoms in the sample. 


As it was shown in our previous work 2^ due to stochastic nature of the electronic 
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damage process the expression ([T]) can be written as a sum of two terms 

(-^(q)) =-^w(q) +-fB(g), (2) 

where 

^w(q) = ^0 exp[fq(Rj - Ri)] (3) 

is the coherent signal containing the structural information and 

= (4) 

i 

is an incoherent background without structural information. Here Iq = J J{t)dt is the total 
fluence of the x-ray pulse. The matrix Wij{q) and vector Bi{q) are defined by the time- 
dependent average values {fi{q, t)) and pulse to pulse fluctuations 6fi{q, t) of the form factor 
fiiQjt) = ifiiQjt)) + of each individual atom through the following relations 

^iM) = Yoj 

Bi{q) = yJ ’ (6) 

where spherically symmetric form factors were considered. The structural term /w(q) © 
determines the degradation of the diffraction pattern due to evolution of the form factors 
while the photoionization process and the background term /b(q') dl]) adds an additional 
background that is due to fluctuations of the individual form factors during the same process. 

An additional background contribution comes from inelastic (Compton) scattering. The 
Compton signal at the detector is given by 

Icoraptoniq) = / J(t) {Si{q, t))dt, (7) 

i 

where {Si{q,t)) is an averaged incoherent scattering function of the atom i (see for details 
Appendix [B]) and brackets have the same meaning as before. Equations ([2]) - ([7]) were used 
in our simulations of diffraction patterns from a biological sample. 

B. Rate equation implementation 

To determine the time-dependent average form factors {fi{q,t)), their fluctuations 
|(5/i(g, and the average incoherent scattering function {Si{q,t)) for each constituent 
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atom type i, a rate equation approach [27| was implemented. First, we define a set of states 
that the atom can potentially occupy. As such states, we consider the electronic ground 
states for all possible occupations of the electronic shells. For example, the carbon atom can 
have between zero and two electrons in each of the Is, 2s, and 2p shells, yielding a total of 
27 states. The time-dependent occupation probabilities p^-i(t) for the ^-th state of the atom 
were obtained by solving a set of coupled differential equations 




( 8 ) 


Here, R^rj-^i{t) denotes the total time-dependent rate of transition from state p to state ^ for 
atom type i. We assume that initially all atoms are in the ground state. The solution of the 
differential equations ([8]) yields the time-dependent occupation probabilities p^;i(t) for each 
state ^ of the specific atom. 

The form factors for each state ^ were obtained from electronic wave functions 

calculated within the Hartree-Fock-Slater (HFS) approximation 28|]. Within this model, 
the average form factors {fi{q,t)) and their fluctuations {\6fi{q,t)\) are given by 

ifiiqR)) = , (9) 

{\suq,t)\^) = {\Mq,t)\^) - l(/.(g,t))r = - !(/.(?,t))r • (10) 

Note that in Eqs. fl9H- ffT0l) we used explicit state-dependent form factors without the addi¬ 


tional assumption that they scale with the number of bound electrons as in Ref. 


29]. In 


the frame of our approach, the valence shells contract significantly on ionization of the core 
electrons, and, as a consequence, the corresponding form factors expand in reciprocal space 
(see Ref. 22]). 

The average time-dependent inelastic scattering function {Si{q,t)) can be calculated in a 


similar way 


^1 


{Si{q,t)) = ^p^;i(t)%(g,t) 

where S^-i{q,t) is the inelastic scattering function for the state ^ 


( 11 ) 




S(A<i) = Z(, - E 


( 12 ) 


r=l 
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Here Z^.^i is the total number of electrons in the atom in state and fl{q) is the form factor 
of the r-th electron in the atom. Here we neglect effects that may forbid excitation of an 
electron from one orbit to another due to Pauli exclusion principle. 

The time-dependent total transition matrix R(f) contains contributions from several elec¬ 
tronic processes 


R(t) = RP^°*°(t) -f + RCompton^^^ ^ + R*’'^P(t) , 


(13) 


where RP’^°*°(t) is the rate of direct photoionization, is the Auger decay rate, R®®'^'^P®(f) 

and R*’'®'P(t) are the rates of secondary ionization produced by escaping and tra ppe d elec¬ 
trons, respectively. These four terms have been considered in our previous work 22j. Here 
we also take two additional ionization channels into account, namely shake-off processes with 
the rate and ionization due to Compton scattering with the rate R‘"°™P*°'^(f), which 

can be important at high x-ray energies. Notice that the first four terms are purely atomic, 
while the latter two are collective effects. See appendix [B] for details on the evaluation of 
the rates. 

In the model of electronic transitions used in this paper we assume that the electron 
plasma thermalizes instantaneously, i.e., the thermalization process is much shorter than 


the pu 
pulses 


se duration. This is considered as a good approximation for comparatively long 


23l |. For x-ray pulses as short as 100 as it is necessary to investigate this question 


in more details. It is well established that non-homogeneous trapped electron gas is formed 


on very early stages of x-ray pulse particle interaction 


23 


3l[ |. To estimate these relaxation 


times the following arguments are typically used. Electrons emitted from a center of the 
spherically symmetric particle of radius R are trapped if their kinetic energy is lower than 
Efrap = e^R'^n/3eo, where e is the electron charge, n is the charge number density and eo is 
the permittivity of vacuum. For Auger energies about ~ 250 eV for the particle with 

o _ 0 

the radius of i? = 15 nm we obtain the charge number density n 2 ■ 10 A at which 
Auger electrons are trapped by the ionized particle. Such small charge density corresponds 
roughly to 3 ■ 10“^ electrons being ionized per atom. Assuming that the dominant ionization 
process is the direct photoionization for the flat-top x-ray pulse we get an estimate for the 
charge number density n{t) ~ a'P^°^°natF{t—to) jT, where Uat is atomic density, F is the pulse 
fluence, and T are times of the pulse start and duration, respectively. From this relation 
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we obtain that the trapping time scales inversely with the photoelectron cross section and 
fluence and is proportional to pulse duration. Our estimates show that for all pulses below 
1 fs and fluences considered in this paper the formation times of non-st at ionary trapped 
electron plasma are below 10 as. 

At the same time, thermalization process of this non-homogeneous, trapped electron gas 
takes place on much longer time scales. Calculations of characteristic thermalization times 


performed according to 


32| give an estimate of about few femtoseconds. By extending our 


model to shorter pulse durations we slightly overestimate the ionization rate of trapped 
electrons. However, this contribution at pulse durations below 1 fs is already significantly 
lower than the contribution from other ionization processes. This is due to the fact that for 
very short pulses Auger electrons do not contribute to ionization process, while secondary 
ionization by escaping and shake-off electrons are producing only low energy secondary 
electrons that can not effectively ionize. By these arguments we can extend our model to 
times as short as 100 as, keeping in mind that we still slightly overestimate the contribution 
from secondary ionization. 

In our model, we assume that the lowest non-vanishing order perturbation theory (LOPT) 
is valid for high energy x-rays in the range of powers up to 10^® W/cm^ and for pulse 
durations down to 100 as. This assumption is based on the fact that ionization for high 
photon energies is well described in the frame of LOPT, if the pulse duration is significantly 
larger than the field period (see e. g. Ref. (33|). 


III. RESULTS AND DISCUSSION 


To analyze the effect of electronic radiation damage, we simulated SPI experiment as 
sketched in Fig. [H For the sample, we used a human adenovirus penton base chimera shell 
34l |. It has a dodecahedral shape with a diameter of 27 nm and contains about 200 000 
nonhydrogen atoms, giving a mass density of about 0.5 g/cm^. To account for typical virus 
densities we have increased this mass density value by three in ionization calculations. The 
ratio between carbon, oxygen and nitrogen atoms in this sample is approximately 3:1:1. 
In our simulations, we neglected the contribution from hydrogen and sulfur atoms to the 
ionization dynamics and the scattering. 

We performed simulations for photon energies of 3.1 keV and 12.4 keV with the exper- 





imental parameters listed in Table [H The lower energy is experimentally attractive due to 
higher elastic scattering power, though at the cost of a lower resolution (about 10 A), while 
the higher energy would be required for reaching a few Angstrom resolution. 

It is important to note that for our model and sample, the photon energy does not affect 
the qualitative ionization dynamics. Increasing the photon energy has two major effects: the 
cross section of the photoionization is rapidly decreasing, and escaping photoelectrons have 
a higher kinetic energy, thus ionizing fewer atoms on their way out. The latter ionization 
process is not dominant for our sample, and the former process merely leads to an effective 
rescaling of the fluence. For light atoms (C, N, O) considered in this work the photoionization 
cross section at 3.1 keV is two orders of magnitude larger than at 12.4 keV. Consequently, the 
ionization dynamics at the photon energy of 3.1 keV and fluence F are similar to dynamics 
observed at 12.4 keV and a fluence 10^ • F. 

All diffraction patterns were simulated in kinematic approximation using mode decom¬ 
position described in Ref. [22 1. 


A. Ionization dynamics 

To describe the contribution of a specihc ionization process a, we hrst introduce sev¬ 
eral quantitative measures. We dehne the probability R'^{t)dt that an atom of the type i, 
undergoes a state change due to a process a during a time interval [t, t -|- dt] as 

= ( 14 ) 

Since all state changes in our model lead to the removal of a single electron, Rf{t) is also 
the rate of ionization. 

To get a global measure of the ionization process, we can integrate these ionization rates 
over time to get the number of ionized electrons due to process a, 5^{t) = Rf{T)dT. 
Normalizing this quantity to the number of electrons of the neutral atom, Zi, and integrating 
with the normalized pulse shape J(t), gives the average degree of ionization. 



The average degree of ionization is a quantity between 0 and 1 that determines how many 
electrons are lost due to process a. The weighting with the pulse shape guarantees that 
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ionization is counted only during the time of the pulse propagation through the sample. For 
example, if Auger ionization starts after the pulse has propagated through the sample, it will 
yield insignihcant contribution to Af, though the degree of ionization can be substantial. 
For our analysis we also introduced quantities averaged over all atoms 


C,N,0 

R“(t) = 5^ 

(16) 

2 

C,N,0 

A“ = ^ M;.A“, 

2 

(17) 

nft = 1 - ^ A“ 

(18) 


a 


with the weights wc,w-^,wo of 3/5, 1/5 and 1/5, corresponding to contribution of C, N, 
and O, respectively. The quantity fib can be interpreted as the normalized average number 
of electrons bound to an atom during the pulse propagation, with rTf, = 1 denoting an 
undamaged atom. 

The time-dependent rates {t) for different ionization processes and average degrees of 
ionization are presented in Fig. [2] and Fig. |3l respectively. In both cases, results are 
shown for the two photon energies of 3.1 keV and 12.4 keV and the same fluence of 10^^ 
photons/pm^. Since the photoionization cross section drops by two orders of magnitude at 
the higher photon energy, this fluence corresponds to two different ionization regimes. At 
3.1 keV, most of electrons are removed from their atoms, so we call this the strong ionization 
regime. In contrast, at 12.4 keV about half of the electrons remain bound even at the end 
of the longest pulse, therefore we call this the weak ionization regime. 

From the data presented in Fig. [2] and Fig. [3] we can conclude that: 

1. Only photoionization and impact ionization by trapped and escaping electrons con- 
tibute substantially to the direct ionization of the atoms. 

2. The net effect of photoionization is independent of the pulse duration (Fig. [3]). Since 
the rate is a time derivative of the number of ionized electrons, the photoionization 
rate is inversely proportional to the pulse duration (Fig. [2]). 

3. As expected, for short pulses of 0.1 fs Auger process is reduced by two orders of 
magnitude in comparison to 10 fs pulses (see Fig. |3]). 
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4. Impact ionization from trapped electrons is the dominant ionization process at rather 
long pulse durations of 10 fs. If the pulse duration is reduced, this process is suppressed 
for the strong ionization regime (Fig. |21(a,b)), and delayed to the end of the pulse for 
the weak ionization regime (Fig. [21(d)). In both cases, the average degree of ionization 
from trapped electron ionization decreases with decreasing pulse duration (Fig. [3]). 
While this decrease is particularly signihcant for sub-fs pulses, we point out that the 
total degree of ionization (see black curve in Fig. [3]) already goes down by half if we 
reduce the pulse duration from 10 fs to 1 fs. 

5. With decreasing pulse duration, impact ionization by escaping photoelectrons becomes 
a relevant process (Fig. [3l [2]). This process becomes especially important for sub-fs 
pulses, and it substitutes the ionization from trapped electrons in the strong ionization 
regime (Fig. [2](a,b)). 


These hndings can be explained from basic considerations. Photoionization in the x-ray 
energy range is for all practical purposes an instantaneous process that only depends on 
the pulse fluence. Hence, any reduction of the pulse length with a constant fluence leads 
to a corresponding increase in the photoionization rate without changing the ionization 
dynamics. Note that the cross section for photoionization of valence electrons is an order of 
magnitude smaller than that of core electrons. This effect causes the apparent shift between 
the photoionization rate and the pulse shape in Fig. [2l(a-c), where all core electrons are 
ionized at the onset of the pulse. While this difference in cross sections can in principle be 
used to create hollow atoms [271], this does not play a role at these particular fluences. 

The impact ionization by trapped electrons is hindered by three factors: First, the cross 
section decreases by about a factor of three for each additional charge of the atom, hence, 
impact ionization becomes a slow process for highly charged atoms. Second, as the atoms 
are ionized, the binding energy of the valence electrons increases rapidly. Finally, the energy 
of the trapped electrons is replenished only by the Auger process. These have a typical 
lifetime of several fs that increases further if there are fewer valence electrons to hll the core 
holes. A decrease of the pulse duration therefore allows to outrun the Auger decay, which 
makes the impact ionization by trapped electrons as negligible. 

We note also that the double-peak form of the trapped electron ionization rate in 
Fig. [2Kc,f) arises from an interplay of these factors. The ionization rate drops initially 
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because the trapped electrons cannot supply sufficient energy to ionize further atoms. At 
that point, Auger decay sets in, leading to a second maximum, after which the sample 
becomes so strongly ionized that the impact ionization becomes inefficient. 

For our sample consisting of light atoms, and for the considered photon energies, ion¬ 
ization from escaping photoelectrons is not as efficient, because the impact ionization cross 
sections drop rapidly with increasing electron kinetic energy. At most about every sec¬ 
ond photoelectron ionizes an atom on its way out of the sample. For long pulses of 10 fs 
(Fig. [2](c,f)), the atoms have already been strongly ionized by the trapped electrons when 
the maximum of photoelectrons are produced. However, for sub-fs pulses, the hrst photo¬ 
electrons encounter a sample of neutral atoms, making subsequent impact ionization more 
likely. Hence, this process plays a role only for the shortest pulse durations 0.1 fs. In the 
strong ionization regime it also appears only at the onset of the pulse (Fig|2](a)). 

As a rough measure of the resulting scattering power of the sample, we can consider the 
square of the average number of bound electrons, n^, Eq. (1T8|) . The results for both photon 
energies and different pulse durations and fluences are shown in Fig. 01 The horizontal dotted 
line shows a cutoff where the sample retains approximately 10% of its scattering power. For 
3.1 keV (Fig. Ufa)) this cutoff is crossed at fluences from 10^^ to 10^'^ photons/pm^ for pulse 
durations from 10 fs to 0.1 fs, respectively. The same behavior is observed for 12.4 keV 
photon energy (Fig. 01(b)) where the cutoff is reached for fluences in the range from 10^^ to 
10^® photons//im^. Altogether, we hnd that reducing the pulse duration from 10 fs to 1 fs 
signihcantly reduces the electronic radiation damage. A further reduction to 0.1 fs yields 
another, but noticeably smaller reduction. 

B. Elastic and inelastic scattering 

As discussed in Section III A1 the hnal signal at the detector has three contributions: 
the elastically scattered coherent signal /w(h)) incoherent background Ib{q), and Comp¬ 
ton background /compton(<?)- Only the coherent signal /w(h) carries information about the 
particle internal structure. It is therefore important to understand how electronic damage 
influences the coherent signal and background contributions. 

The coherent signal as a function of the momentum transfer for different fluences and 
pulse durations is shown in Fig. [5l The intensity was calculated according to Eq. ([3]) and 
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angularly averaged over all detector pixels of constant |q|, giving the angular averaged 
intensity, (Jw(q))<^ = (27r)“^ Iwil, per Shannon angle. 

If we disregard the background contribution , we can dehne the maximum achievable res¬ 
olution by requiring a minimum of 10“^ photons per Shannon angle for successful orientation 
{q, 10|. In practice, this number may be higher due to artifacts and noise. At 3.1 keV photon 


energy and a fluence of 10^^ photons/pm^, a pulse duration of 10 fs allows to achieve 8 A 
resolution, while 1 fs and 0.1 fs pulses allow to reach 4 A. For 12.4 keV and a fluence of 10^^ 
photons/pm^, we can achieve about 3 A resolution; a further increase towards 1 A is only 
possible by increasing the fluence even further and at the same time having pulses of less 
then 1 fs duration. Analysis of the results presented in Fig. [S] shows a substantial difference 
between the strong ionization regime Fig. O^a-c) and weak ionization regime Fig. El^d-f). In 
the former case scattered intensities are substationally lower than intensities corresponding 
to an undamaged sample even at very short pulse durations of 0.1 fs (see Fig. [5](a)) and in 
the latter case they are very close to each other. Another important effect is the saturation 
of the scattered intensity. At high fluences in the strong ionization regime an increase of 
the XFEL intensity by one order of magnitude does not lead to the same increase of the 
scattered intensity. 

In the following, we consider a resolution of 10 A for 3.1 keV and 3 A for 12.4 keV, 
corresponding to go ~ 0.6 A“^ and go ~ 2 A“^, respectively. Fig. [6] shows the averaged 
scattering intensities {Iw{Qo))cI> at these momentum transfer values as a function of pulse 
duration and fluence. 

For an undamaged sample. Fig. |6] shows a linear relationship between the incoming pulse 
fluence and scattered intensity. However, due to ionization of the sample, this relationship 
breaks down at high XFEL intensities. For 10 fs pulses, deviation from the linear scaling law 
starts at 10^^ (10^“^) photons/pm^ for 3.1 (12.4) keV, and becomes signihcant at one order of 
magnitude higher fluence. This is caused by a substantial decrease of the number of bound 
electrons that can scatter (see Fig. 0]). Reducing the pulse duration to 1 fs increases the 
scattered intensity by reducing electronic radiation damage. A further reduction to 0.1 fs 
gives another, but considerably smaller increase, in qualitative agreement with the radiation 
damage observables (Fig. 0]). 

The results here put both lower and upper boundaries on acceptable fluences for imaging 
the test particle. To achieve the required resolutions, a minimum fluence of 2 ■ 10^^ (10^'^) 
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photons/yum^ for 3.1 (12.4) keV is strictly required to get enough scattered signal. At 10 
fs pulse duration, however, an increase in fluence no longer translates into an increase in 
scattered intensity for intensities beyond 10^^ (10^^) photons/pm^. An increase of the fluence 
by one order of magnitude increases the scattered signal only by a factor of two. Decreasing 
the pulse duration to 1 fs already improves the scaling significantly, with another smaller 
gain when going to 0.1 fs pulses. Still, even for the shortest 0.1 fs pulses, there is little 
advantage from increasing the fluence beyond 10^^ (few 10^®) photons/pm^ for 3.1 (12.4) 
keV. 


As a simple measure of the background effects, we can use the ratio between the respective 
background contribution (incoherent signal /b(q') or Compton background /compton(Q')) and 
the coherent signal. 


rB/Compton(<?) 


/ Compton ( O') 

(/w(q))^. 


(19) 


This measure for the incoherent background Tsiq) is shown in Fig. [71(a,c) for both photon 
energies and 1 fs pulse duration. Note that the oscillations in these hgures are caused by 
the speckle pattern of the diffraction image; the background /b(q') is a smooth function 
of q. We considered a cutoff of 10% shown as a horizontal dashed line, after which the 
background becomes a signihcant feature of the diffraction pattern and complicates the 
analysis, especially the orientation of the single-shot diffraction patterns. 

The dependence of rB((3'o) on the XFEL fluence and pulse duration is presented in 
Fig. [71(b,d). The background rises continuously from negligible noise to the dominant con¬ 
tribution as the fluence increases. For the highest fluences, FB(go) shows saturation for all 
x-ray parameters, at values up to one. A reduction of the pulse duration slightly reduces the 
background and the saturation value of FB(go)- For 10 fs pulses, the cutoff is reached for 
fluences of 3 ■ 10^^ (10^®) photons/pm^ for 3.1 (12.4) keV x-rays. Reducing the pulse dura¬ 
tion further to 0.1 fs increases the allowed fluences to about 2 ■ 10^^ (6 • 10^®) photons/pm^. 
Hence, the use of very short pulse durations is experimentally attractive to suppress this 
background contribution. Note that the restrictions on the pulse fluence are similar to those 
from considering only the coherent scattering /■w(q). 


A comparison of the coherent signal and Compton scattering contribution is shown in 
Fig. [8] for pulses with 1 fs duration, different fluences and both photon energies. Only the 
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contribution from bound electrons was taken into account, the contribution from trapped 
and escaping electrons was neglected. Hence, the presented results could be considered as 
a lower boundary. As expected the Compton scattering becomes more important at high 
photon energy, and dominates the signal for high momentum transfers at g > 3 A“^. The 
Compton scattering is relatively weak for soft x-rays, never reaching the coherent signal. 

The ratio Tcompton (Eq. flT^ l for different pulse parameters is shown in Fig. |9l We point 
out that in our simulations, the relative Compton background from bound electrons is larger 
than in 2^. This increase is due to two factors. The atoms in our simulations are stronger 


ionised on average due to inclusion of the electron impact ionization that leads to a stronger 
suppression of the coherent scattering. Also, the explicit inclusion of nitrogen/oxygen atoms 
with more valence electrons (3 and 4 respectively versus 2 for carbon) increases the Compton 
scattering contribution, since this process dominantly occurs on weakly-bound electrons. 

We found that the limit of 10% background is always surpassed for the momentum transfer 
values larger than g > 1 (g > 1.5 A“^) and photon energies of 3.1 keV (12.4 keV) (see 
Fig. ini(a,c)). We also observed that the dependence of the Compton background at a constant 
momentum transfer value go on the pulse parameters is rather weak (see Fig. [ni(b,d)). The 
Compton contribution practically does not depend on the pulse duration and Fcompton(go) 
increases only by a factor of two to three for fluences above F > 10^^(10^®) photons/pm^ at 
3.1 (12.4) keV photon energy, respectively. Effectively, considering the maximum acceptable 
background level to be 10%, the Compton scattering limits the achievable resolution to 
approximately 6 A (4 A) for the soft (hard) x-rays. 

Our simulations show that the Compton scattering gives a substantial contribution in 
the hard x-ray scattering conditions and less important in the soft x-ray range. Without 
a proper treatment of this background or use of energy-resolved detectors, few Angstrom 
resolution limit will be difficult to reach for small non-crystalline particles. 


IV. SUMMARY AND CONCLUSIONS 


In summary, we have extended our previous approach 


22| on ionization dynamics of 


biological samples to incorporate shake-off ionization and Compton scattering. We studied 
the ionization dynamics and the scattered signal for ultrashort XFEL pulses from 0.1 fs to 
10 fs. We used an adenovirus shell as a test sample, and considered soft(hard) x-ray pulses 
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with 3.1(12.4) keV photon energy and 10(3) A target resolution, respectively. 

By introducing appropriate measures, we quantihed the contribution of the single ion¬ 
ization mechanisms to the electronic radiation damage. In particular, we found that with 
sufficiently short pulses (on the order of 1 fs and less) it is possible to outrun the ionization 
from the trapped electron gas and therefore to reduce the electronic damage significantly. 
Our simulations show that it is not necessary to use sub-fs pulses; a considerable damage 
reduction is already realized for pulse durations of 1 fs. 

To translate this into a useful fluence limit, we also analyzed the scattered intensity, 
which has three contributions. These are a coherent signal that contains all the structural 
information, an incoherent background that is due to statistical fluctuations of the form 
factors of individual atoms, and the Compton (inelastic) background. For a given sample, 
they put different boundaries on the XFEL fluence and achievable resolution in different 
ways. If the XFEL fluence is too large, most of electrons are striped from the atoms and the 
sample does not scatter anymore. In this case scenario the coherent signal does not increase 
any more with the increasing XFEL pulse fluence. 

We found that for a given XFEL fluence the use of shorter pulses increases the coher¬ 
ent signal and reduces the incoherent background. In concord with the radiation damage 
measures, we found a large improvement already by reducing the pulse duration from 10 fs 
to 1 fs. Even with the shortest pulses of 0.1 fs, however, there is an upper limit of about 
10^“^ (few 10^^) photons/pm^ for 3.1 (12.4) keV photon energy. Also at these fluences the 
incoherent background becomes relevant, the signal-to-noise ratio drops down, and recovery 
of a structural information will become a challenging task. 

We also studied the Compton background, and found that it provides a strong contribu¬ 
tion to the scattering signal at a few Angstrom resolution. For our test sample, we found 
that Compton scattering limits the available resolution to about 6 A (4 A) for a photon 
energy of 3.1 (12.4) keV. This background is virtually independent of the pulse parame¬ 
ters, and cannot be suppressed by the use of short pulses or small fluences. To overcome 
this problem a special study of the Compton contribution and, possibly, dedicated energy 
resolution detectors will be necessary in the future. 

Alltogether, our analysis shows that SPI experiments are still challenging, especially for 
small biological samples of 30 nm size and below. To reach subnanometer resolution we 
suggest to use pulses of about 1 fs and the XFEL fluence that is below high ionization 
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regime. In order to obtain access to high resolution a substantial amount of diffraction 
patterns should be accumulated at these XFEL conditions. Our results show that it is 
not advisable to go to high ionization regime of XFEL operation since ionization dynamics 
prevents substantial increase of the scattered intensity with the raise of the XFEL power. 
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Appendix A: Elastic scattering of extremely short pnlses 

Here we present the general formalism of elastic scattering of ultrashort FEL pulses on a 
single particle in kinematical approximation. 

We consider the amplitude of the incoming x-ray wave in the form 


Ein{r, t) = Ein{r, t) exp(iki ■ r - iut). (Al) 

Here E’j„(r, t) is the slowly varying amplitude, u and ki = (27r/A)ni are the average frequency 
and momentum of the incoming x-ray wave with the average wavelength A. The incident 
direction of the wave is defined by the vector ni. 

In the frame of the hrst Born approximation the instantaneous amplitude of the outgoing 
wave at the time t at the detector position coordinate u for a narrow bandwidth light can 
be presented as |35j 


Eout{u,t) = ^ f drp(r,t - (A2) 

%X J ^ur 

where = Rur/c is the time delay for the light propagating from the position r in the 
sample to the position u in the detector and R^r is the distance between these two points 
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(see Fig. Al). Substituting now expression flAlj) into Eq. flA2j) we find for the amplitude of 
the scattered wave 


■^out 


(U,t) = \ [ drp(r, t - Tr) g.k,r-.^(t-r.)_ 


i\ 


K, 


(A3) 


In the Fresnel limit we have for the distance Rur between the points r and u 


Rur = L — Ilf ■ T + 


U 


2L 


(A4) 


where L is a distance between the sample and the detector and vector nf = kf/|k| defines 
the direction of the outgoing wave. Assuming far-field limit {(P/\L « 1, where d is the 
size of the sample), as is typical in the case of the single particle imaging experiments and 
substituting flA4p in Eq. flA3p we find for the scattered field 

pikL—iujt 

EoutiuR) = 


A{q,t) = / drp{r,t - Tr)Ein{r,t-Tr)e 


(AS) 


where the momentum transfer vector q = (kf — ki)+qu with qu = k{u/L) is introduced. This 
is the far-field expression for the instantaneous value of the scattered field at the detector 
position u. 

We write now the electron density of a single particle in a usual way as a sum of instan¬ 
taneous electron densities Pi{r,t) of each atom at the position Rj 

N 




(A6) 


2=1 


We assume now that for the sufficiently short femtoseconds pulses considered here atomic 
positions do not change during the pulse propagation and all time dependencies are due to 
electronic changes in the individual atoms consisting the particle. We also neglect here all 
cooperative effects. Substituting expression flA6p for the electron density into Eq. flA5D a nd 
performing the change of variables we obtain for the amplitude of the scattered field 441 


N 


A(q, t) = / dr'pi{r', t - TR,)Ein(Ri + r',t- TR^)e (A7) 

i=i 

Assuming here that the incident x-ray field is uniform on the size of a single atom 
Ein(Ri + r',t - Tr^) ^ Ein{Ri,t — tr^) and introducing the time-dependent atomic form 
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factors /i(q, t) = / pi{Y,t)exp{—i(\ ■ r)(ir we obtain the following general expression for the 
scattered amplitude flA7D 


N 

dl(q, t) = ^ “ rR.)Eir,{Ki, t - (A8) 

i=l 

This expression for the scattered amplitude from a single particle differs from the tradi¬ 
tionally used kinematical expression by two important features. First, it contains the time 
delayed atomic form factors /i(q, t and, second, it contains the instantaneous time 

delayed incident held Ein(Ri,t — t/jJ at each atomic position. 

For a typical single particle imaging experiment diffraction pattern recorded on the de¬ 
tector will be given by the intensity of the waveheld dehned by Eq. flASp . As soon as present 
detectors do not have femtosecond time resolution the measured signal will be, necessarily, 
integrated over the time of the pulse duration and will be given by 

N 

J(q) = / /(q,t)dt = ^ e-*q (R.-Rdx 

i,k=l 

X J dtf*{q,t - TR^)fkiq,t - TR^)E*^{'Ri,t - TR^)Ein{Rk,t - Tr^). (A9) 

As soon as this experiment will be repeated many times we would have to average results 
of these measurements. This is equivalent to ensemble averaging of time-integrated intensity 
distribution flA9p 

r ^ 

< /(q) > = / < /(q,t) > dt = ^ g-*q-(Rfc-Ri)^ 
d i,k=i 

X J dt< /;(q, t - TR.)fk{q, t - rRjE;„(Ri, t - r/jjEi„(Rfc, t - Tr^) > . (AlO) 

Assuming now that huctuations of the incoming waveheld are statistically independent 
from the huctuations of the electronic system we can factorize averaging in Eq. flAlOp into 
a product of two terms 


N 


< J(q) > = ^ e 

i^k=l 


X 


X j dt < TR.)fkiq, t-TR^)> X 

X < E*^{Ri, t - TR,)Ein{Rk, t-TR^)> . 


(All) 
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At that stage we introduce for the incoming x-ray held a mutual coherence function 361 
rin(ri, r2; h, ^ 2 ) = t - ti)Ein{r2; t - ^ 2 )) (A12) 

and its normalized version the so-called complex degree of coherence 


7m(ri, r2; ^ 1 ,^ 2 ) = ri„(ri, r2; ^ 1 ,^ 2 )/ J{r2;t), 
where the incident intensity is dehned as 

J(r;t) = ri„(ri = r2 = r;ti = ^2) = (|-Ei„(r,t)|^) . 
Substituting these dehnitions in Eq. flAllD we obtain 

N 

< /(q) > = ^ g-*q'(Rfc-Ri); 


(A13) 


(A14) 


X 


i,k=l 


X J dt < rHj/fc(q, t - tr^) > X 

where nk = tr, - tr^. 

This is very general expression for the averaged intensity in a single particle imaging 
experiments derived in kinematical approximation. It takes into account degradation of 
contrast of the coherently scattered intensity from a single particle due to two effects. First, 
due to time evolution of the electronic structure of each atom as a result of fast ionization 
while propagation of the femtosecond x-ray pulses. Second, due to partial coherence (spatial 
and temporal) of the incoming radiation. As it was demonstrated in a series of experiments 
371 . l38| x-ray pulses from XFEL sources have a high degree of spatial coherence and limited 
temporal coherence that could, in principle, degrade the contrast of the scattered intensity 
at high resolution. 

At the next stage we will assume that the incoming radiation is fully coherent spatially and 
temporally, we will neglect retardation effects, and consider that the intensity distribution 
is spatially uniform over the size of the sample. This leads to the following expression for 
the scattered intensity 


N 


< J(q) > = ^ e R-i) 


X 


i,k=l 


X / dtJ{t) < /;(q,t)/fc(q,t) > , 


(A16) 
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This expression was used as a starting point for the analysis in our previous paper [22| and 
in the present work. 


Appendix B: Details on the rate calculation 

The direct photoionization rate is calculated as j (t), where is 

the photoionization cross section 45| . Since in the considered energy range the cross section 
is a slowly varying function of energy we use the cross section at the central frequency of 
the x-ray pulse. The photoionization cross sections as well as the Auger rates 

were calculated within the HFS approximation 2^; the explicit expressions can be found, 
for example, in [^ . 

During photoionization of core electrons, especially at large photon energies, the screening 
potential which is felt by other electrons changes abruptly which can lead to the further 


emission of one or more electrons in a so-called shake-off process 


39[ |. If the photoionization 


transition occurs from the state rj to the state (, it can be accompanied by the shake-off 
transition to the state ^ with the probability Then the rate of such transition is 


p|^shake^photo(^) total rate of various shake-off transitions to the same state ^ is 




shake ^photo 


C C»7 Cv 

The Compton ionization rates are defined as J(^ Th e Compton 

cross sections 5 -Compton calculated in impulse approximation as 

duKN{d) 


30 


m 


a. 


Compton _ 


-Sm-i{q)dn . 


(Bl) 


J dVt 

Here daKN{0)/dfl is the Klein-Nishina differential cross section, 9 is the scattering angle, 
and is the incoherent scattering function for electrons at the m-th shell for an atom 

in initial state The cross section describes the Compton process in which an 

electron is ejected from the m-th shell of an atom in the state For low photon energies 
(e = E^/nieC^ 1) the Klein-Nishina differential cross section within a good approximation 


can be simplified to 


42| 


da kn{9) 
dVt 


rl[l -F e(l - cos 6 ')] P{6) , 


(B2) 


where rg is the classical electron radius and P{9) is the polarization coefficient. In our 
simulations for photon energies up to 12 keV we used this approximation for Klein-Nishina 
cross section. 
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The incoherent scattering function for electrons at the m-th shell Sm-,^{q) was calculated 


in the frame of HFS approximation as 


30| 






(B3) 


2 = 1 


where f^{q) is the form factor of Tth electron on the m-th shell, and the number of 
m-shell electrons in state Here we neglect effects that may forbid excitation of an electron 
from one orbit to another due to Pauli exclusion principle. 

It is well established (see for example Refs. 2^, l23j) that secondary ionization can 


significantly change the ionization behavior in the biological particle. In our model we treat 


the contribution of the secondary ionization similar to our previous work 


22 ] (see also 231) 


and neglect the details of the trapping process as a function of the particle charge and 
electron dynamics. However, our model reproduce well most of the important features of 
the secondary ionization discussed in Ref. 231. 

We assume that photoelectrons originate in the center of the biological particle and during 
their escape produce secondary electrons through impact ionization. The biological particle 
is considered to be in the form of homogeneous sphere with the radius R and volume V. 
The rate of secondary ionization produced by the escaping high-energy photoelectrons is 
calculated according to j22l. l23| 

R AM , . 

(B4) 


^escape _ ^impact ^-^photo 


where diVphoto/dt is the production rate of photoelectrons. The impact ionization cross 
sections were calculated using the binary-encounter Bethe model 43| with electron 

orbital parameters obtained from the HFS calculations. When a secondary electron is pro¬ 
duced, its kinetic energy was set to a constant value Eq = 25 eV [23 1. 

In our model Auger, shake-off, Compton as well as all secondary electrons produced by 
electron impact ionization were supposed to be trapped by the positively charged particle 
and thermalized instantaneously into a homogeneous electron gas with a Maxwell-Boltzmann 
distribution of velocities. The temperature T of the gas was determined by the average 
kinetic energy of the trapped electrons as E^^ap = (3/2)fcT, where fc is a Boltzmann constant. 
The rate of collisional ionization produced by these trapped electrons was determined by 

Rt-P(t) = ntrap(t) , 
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where we specifically introduced dependence on electron velocity v in impact ionization cross 
section Here, ntrap(^) is the time-dependent density of trapped electrons and Ot 

is thermodynamic average with the appropriate Boltzmann factors. 

The time evolution of the trapped electron density ntrap(^) and total kinetic energy ii^trap(^) 
of the trapped electron gas were calculated from the set of equations 
dut 


^trap 

dt 


i Cr; 


E E <T' + «""'(*) + (‘) + )PnMe., (B5) 


dE, 


trap 


dt 


i 


E E (O-Bg: 


shake 




(B6) 


Here the index i denotes the atom type (carbon, nitrogen etc.) with the corresponding 
atom density p*. The energy is the energy of the released Auger electrons, which 

is assumed to be approximately a constant ^ ^Auger fQ^- particular 

atom. The energy of the released shake-off electrons is assumed to be about 1.8 times 

the binding energy of the electron 39|]. The energy of the electrons ionized by Compton 
scattering ^Compton assumed to be constant for a particular atom type. The 

energy is the binding energy of the electron that is released through impact ionization, 
and Eq = 25 eV is the energy assigned to secondary electrons produced by photoelectrons. 
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Figures 



FIG. 1: (Color online) Schematic of a single particle coherent X-ray diffractive imaging experiment. 
A single FEL pulse illuminates the sample of size d from the left and scatters from it, with the 
diffraction pattern measured by a detector of size at a distance L from the sample. 
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FIG. 2: (Color online) Average ionization rates fC{ t) of different radiation processes for photon 
energies of 3.1 keV (a-c) and 12.4 keV (d-f). The pulse durations are 0.1 fs (a,d), 1 fs (b,e), and 
10 fs (c,f). The pulse fluence is 10^^ ph//im^ in all cases. The black thin solid line shows the pulse 
shape. 
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FIG. 3: (Color online) Average degree of ionization A for different ionization processes and 
pulse durations for photon energies of 3.1 keV (a) and 12.4 keV (b). The pulse fluence is 10^^ 
photons//zm^ in both cases. 




















FIG. 4: (Color online) Square of the average normalized number of bound electrons, for different 
pulse durations and fluences as a figure of merit of the samples scattering power. Photon energies 
are 3.1 keV (a) and 12.4 keV (b). The horizontal black dashed line corresponds to a cutoff of 10% 
of the scattering power of the undamaged sample. 
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FIG. 5: (Color online) Angular averaged coherent signal (/w(q ))0 for photon energies of 3.1 keV 
(a-c) and 12.4 keV (d-f) and different fluences, defined in ph/fim'^ in the insets. Pulse durations 
are 0.1 fs (a,d), 1 fs (b,e), and 10 fs (c,f). The horizontal dashed black line corresponds to the 
requirement of 10“^ photons per resolution element. The curves have been smoothed with a 
Gaussian filter to remove high-frequency oscillations from the speckle pattern. Thin dashed lines 
correspond to undamaged sample and fluence 10 ^^ ph/pm?. 
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FIG. 6: (Color online) Angular averaged coherent signal < Iw{qo) >0 for different pulse durations 
and fluences, defined in ph/fim? in the insets. The dashed black line represents the lower limit of 
10“^ photons per resolution element. The solid black line corresponds to undamaged sample. The 
quantities are evaluated at 10 A resolution for 3.1 keV photon energy (a) and 3 A resolution for 
12.4 keV (b). 
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FIG. 7: (Color online) (a,c) Relative contribution of the incoherent background, rB(g), as a func¬ 
tion of momentum transfer for different fluences, defined in ph /in the insets, photon energies 
of 3.1 keV (a) and 12.4 keV (c), and pulse duration of 1 fs. (b,d) rB((7o) as a function of pulse 
fluence for different pulse duration at a resolution of 10 A for 3.1 keV (b) and 3 A for 12.4 keV 
(d). The horizontal dashed black lines represent an upper acceptable limit of 10% background. 
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FIG. 8: (Color online) Comparison of the angular averaged coherent signal {Iw{fi))(t> (solid lines) 
and Compton signal /compton(<?) (dashed lines) for photon energies of 3.1 keV (a) and 12.4 keV (b) 
and different fluences, defined in ph/nm^ in the insets. The pulse duration is 1 fs in both cases. 
The horizontal dashed black line corresponds to the requirement of 10~^ photons per resolution 
element. 
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FIG. 9: (Color online) (a,c) Relative contribution of the Compton scattering, rcompton(Q') as a 
function of momentum transfer for different fluences, defined in ph/fim^ in the insets, photon 
energy of 3.1 keV (a) and 12.4 keV (c), and pulse duration of 1 fs. (b,d) rcompton(Q'o) as a function 
of pulse fluence for different pulse durations at 10 A resolution for 3.1 keV (b) and 3 A resolution 
for 12.4 keV (d). The horizontal dashed black line represents an upper acceptable limit of 10%. 
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Tables 


TABLE I: Parameters 

used in t 

le simulations. 

photon energy (keV) 

3.1 

12.4 

pulse duration (fs) 

0.1 - 10 

0.1 - 10 

sample-detector distance (mm) 

100 

100 

detector size (mm) 

200 

400 

number of pixels 

320 X 320 

1536 X 1536 

pixels per speckle 

2.4 

1.4 

Achievable resolution at the detector edge (A) 

5.2 

0.95 A 
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